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Abstract 

A phenomenological model of the Landau type forms the basis for a study of 
elastic distortions near the spin-Peierls transition T c in CuGeO^. The atomic 
displacements proposed by Hirota et al. [Phys. Rev. Lett. 73, 736 (1994)] 
are accounted for by the model which includes linear coupling between Cu 
and O distortions. Cu displacements are seen to be responsible for anomalies 
in the elastic properties at T c , whereas incipient O distortions give rise to 
temperature dependence below T c . A discussion of possible critical behavior 
is also made. 
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I. INTRODUCTION 

Due to the nature of its microscopic origins, attention has focused on both magnetic^ 
and structural propertiesUll of the recently discovered spin-Peierls system CuGeO^. It 
is now well established that this material is the first example of an inorganic compound 
which exhibits the simultaineous opening of a gap from a singlet ground state and a lattice 
dimerization in S=l/2 spin chains. In contrast with previously studied organic spin-Peierls 
systems, large single crystals of CuGeOj, can be made which facilitate neutron and X-ray 
diffraction studies in addition to the measurments of crystal-direction-dependent thermody- 
namic properties. Such measurements have revealed large anisotropic strain effects due to 
the dimerization, in contrast with simple quasi-one-dimensional models of Cu displacements. 
Such effects are examined in the present work based on a phenomenological Landau-type 
model which describes the atomic displacements below the transition temperature T c as 
proposed by Hirota et a/.i In contrast with the principal Cu displacements which give rise 
to anomalies at T c , secondary O distortions are found to be important in determining the 
temperature dependence of elastic properties below T c . A linear coupling between the two 
types of ionic displacements is proposed which is electronic in origin. Good agreement with 
available experimental results is demonstrated. This study also serves to compliment and 
extend the brief Landau-type descriptions of some elastic properties given in Refs. |5|,[rT]|,|H)|. 

Of the ten ions per high-temperature (above T c ) orthorhombic unit cell (a = 4.81, b = 
8.47, and c = 2.941 A), at sites R = R/ + Vj (where R; are Bravais lattice vectors and 
j — 1 — 10), there are two Cu and four O ions which displace at T c in the model of Ref || 
These are shown schematically in Fig. 1, with the high-T positional parameters Vj given by 

C Ul : vt = (0,0,0); Cu 2 : v 2 = (0, |, 0) 

O s : v 3 = (x,-y, |); 4 : v A = (-x, y, §) (1) 
5 : v B = (-a:, j - V, §); 6 : v 6 = (x, y + \, |). 

These displacements thus double the unit cell in the a and c directions below T c . A simple 



order-parameter description of the ionic movements below T c may be written as 

p(R) = p j coa(Q-R l ) (2) 

where Q = (|, 1, |) and here j = 1 — 6. The displacements shown in Fig. 1 are then 
described with 

Pi = -PzZ] p 2 = ~P 1 

p 3 = -p x ± - p y y; p 4 = -p 3 (3) 
p 5 = -p x Z + py-y; p 6 = -p 5 

where p x and p y are associated with the O distortions and p z represents Cw-ion displace- 
ments. It is argued below that simple ion-ion interactions derived from symmetry arguments 
can reproduce the configuration of displacements depicted in Fig. 1. A desirable feature 
of any such model is to account for the simultaneous displacement of both Cu and O ions 
(which is also relevant to possible critical behavior). It is for this reason that that ionic 
density-density coupling must also be considered. In addition, simple coupling between p 
and the strain tensor reproduces the principal features of measured elastic properties without 
evoking higher-order interactions as in Refs. ^,|T0|,|T~5 . 



II. LANDAU FREE ENERGY 

Landau-type free energies may be constructed using symmetry arguments^. It is as- 
sumed here that magnetic degrees of freedom are integrated out so that the free energy 
depends explicitly on p only. In the case of orthorhombic symmetry, there are three inde- 
pendent invariants which contribute to free energy density / at second order in the vector 
p, which may be chosen as 

f 2 = D' x p 2 x + D' y pl + D' z p 2 z . (4) 

The independence of these three coefficients requires that p x ^ p y ^ p z , as observed ex- 
perimentally. Alternatively, one can choose the three terms as = Ap 2 + D x p 2 x + D y p 2 , 
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where the A-term represents an isotropic interaction. Minimization of the anisotropy at the 
Cu sites, for example, gives p ± , p 2 || ±z if D x {Cu) > and D y {Cu) > 0. In the following 
analysis, only isotropic terms are considered in detail, which is sufficient for the purposes of 
the present study. 

A. Second-order terms 

A general expression for isotropic two-ion interactions is 

F 2 = W A(R-R')p(R).p(R'). (5) 

Z R,R' 

Using the form (2) for p(R) in this relation yields 

F 2 = ^J2 A AQ)PrP^ (6) 

where 

^■fc(Q) = E A (r + Vj ~ v fc ) cos(Q • r) (7) 
r 

with t = Hi — TL'i. Considered below are the near-neighbor interactions within and between 
the Cu0 2 chains which are likely to be important. We recall that CuGe0 3 is somewhat 
quasi-one-dimensional . 



1. Cu-Cu within chains 

Interactions between near-neighbor Cu ions along the c axis involve j — k — (1,2) with 
An = ^4( c ) cos 7r + A(— c) cos7r, and with a similar expression for A 22 , so that 

F 2A = -A(c)[pi- pi + p 2 ■ p 2 ] 

= -2Aup 2 c (8) 

where An = A(c). 
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2. Cu-Cu between chains 



Here j ^ k = (1, 2) with, e.g., A 12 = A(-v 2 ) + A(-v 2 + by) cos27r, giving 

F 2B = 2A(v 2 ) Pl • p 2 . (9) 

If A(v 2 ) > 0, then this term is minimized by the configuration p 1 || — p 2 , as in Fig. I. It 
can be expected, however, that this term is smaller than F 2 a due to the larger interaction 
distance involved. Its contribution is given by F 2 b = —2A\ 2 p 2 C) where A i2 = A(v 2 ). 

3. 0-0 within chains 

Interactions between the two O ions which surround the Cu of Fig. 1 correspond with 
j ^ k = (3, 4) and (5, 6). Using A 34 = A(v 3 - v 4 ) = A 43 = A Fj6 = A 65 , the energy is 

F 2C = A(v 3 -v 4 )[p 3 - p 4 + p 5 - Pq\. (10) 

If A(v 3 — v 4 ) > 0, then the configuration p 3 \\ — p 4 and p 5 || — p 6 , as in Fig. 1, is stabilized 
by this term, giving F 2C = —2A 3i p 2 0) where A 34 = A(v 3 — v 4 ). 

4- 0-0 between chains 

With j ^ k = (4, 5) and (3, 6), the resulting contribution to the free energy is 

F 2D = A(v 4 -v 5 )[p 4 -p 5 + p 3 -p 6 ]. (11) 

Although parallel or anti-parallel configurations are minimized by this term, it is expected 
to be rather small. Fourth-order terms considered below prefer a configuration p 4 _L p 5 
and p 3 _L p 6 giving F 2D = 0. As emphasized above, however, anisotropy effects will alter 
somewhat this idealized state. 
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5. Cu-0 within chains 



Considered here are the near-neighbor Cu — O interactions, so that j ^ k = (1, 3), (1,4), 
(2,5), as well as (2,6). Due to the fact that neighboring Cu ions along the c axis displace 
in opposite directions, and that the O ions (3,4) and (5,6) are located mid-way between 
these Cu ions, this two-ion interaction term is zero. This is illustrated by considering 
Ai 3 = A(— v 3 ) + A(c — v 3 ) cos 7r = 0, where A(— v 3 ) = A(v 3 ) = A(c — v 3 ) has been used. In 
the absence of Cu — O coupling at second order, there is not a simultaneous displacement 
of both species. This result is not in agreement with the observed behavior. 

6. Cu-0 within chains: density- density interactions 

In addition to the two-ion intereactions which involve the vector displacements 
p(R), there are also ionic charge density couplings which involve the perturbed density 
p(R) = \p(R)\. Such interactions are of the type which have been considered within the 
framework of density functional theoryEI and are of the form 

f g = Wg(R- R')p(R)p(R0 

/ R,R' 

= 2E^*¥* ( 12 ) 

where 

Gji = E^ + vr^)- ( 13 ) 

r 

For simplicity, only the near- neighbor Cu—O coupling is explicitly considered here (the other 
interactions discussed above are simply renormalized by charge density coupling). This gives 
a contribution 

F G = 4:G 13 p c po, (14) 
where (?i 3 = G(v 3 ). Linear coupling of this type yields (see below) po oc po, as expected. 
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Temperature dependence appears through the usual entropy term of mean-field theory^ 

F e = faT5> 2 (R). (15) 

R 

With the simplifying assumption that p x ~ p y (see below), all of the second-order terms 
considered above can be written as 

F 2 = A c p 2 c + A p 2 + Acopcpo, (16) 

where 

Ac = aT-2A n -2A 12 , A Q = aT-2A M , A co = 4G 13 . (17) 

B. Fourth-order terms 

We omit here consideration of fourth-order anisotropy terms and note that, in general, 
isotropic interactions involving the vector p(R) are nonlocal and can be expressed as0 

Fa = \Y. 5(R 1 ,R 2 ;R 3 ,R 4 )[p(i?i) • pCR 2 )][p(i* 3 ) ■ p(# 4 )], (18) 

where S depends only on the relative distances | Ri — Rj | . For simplicity, and for the 
puposes of demonstrating that the present model is plausible, only two types of terms will 
be considered explicitly here: 

4 R 

= B [ P 2 c + 2p 2 o } 2 , (19) 

which comes from entropy considerations,!^ as well as the simple two-ion interactions of the 
form 

F ib = jJ2 B(R - R)[p(R) • P(R')] 2 - (20) 

4 R,R' 

The five near-neighbor couplings as considered above for F 2 yield a result 
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F 4b = B 11 [pi + p*] + B 12 (p 1 -p 2 ) 2 

+ |S 34 [(P3 • P 4 ) 2 + (P 5 ■ P 6 ) 2 ] + |5 45 [(p 4 • p 5 ) 2 + (p 3 • p 6 ) 2 ] 

+ ^i3[(Pi ' Ps) 2 + (Pi • P 4 ) 2 + (P 2 • P 5 ) 2 + (P 2 • Pe) 2 ]- (21) 
Configurations minimized by these terms are as follows: 

For B 12 < 0, p 1 || ±p 2 ; 

/or 5,34 < 0, p 3 || ±p 4 ; p 5 || ±p 6 ; 

/or £45 > 0, p 4 JL p 5 ; p 3 JL p 6 ; 

/or £ 13 > 0, p c JL p D . (22) 

Note that the complete minimization of the term _B 45 occurs only if = p y when (3) is used, 
but that this equality is in competition with anisotropy terms. It appears, however, that 
this condition is nearly satisfied in CuGe0 3 . These results, together with the minimization 
conditions from consideration of F 2) fully explain the displacements of Fig. 1, as described 
by Eqs. (l)-(3). 



C. Equilibrium Properties 
The full free energy is then given by 

F = A c p 2 c + A p 2 o + AcoPcPo + B c p% + B p% + BcopIpI, (23) 

where B c = B + B u + B\ 2 , B = 4_B + -B34, and B C o = 4-B . Due to the nature of the 
spin-Peierls interaction, it is assumed here that Cu and O displacements represent primary 
and secondary order parameters, respectively. Minimization of F with respect to po yields 
a result which can be expressed as an expansion in odd powers of pc, at least near the 
transition temperature, 

Po = ap c + (3p 3 c + ---, (24) 

where 
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a = -\AcojAo- (3 = -a/A (B C0 + \B Q A 2 co / A 2 ). (25) 

This establishes a linear relationship between Cu and O displacements near T c , but that a 
more complicated behavior can be expected at lower temperatures. Using this result in the 
above expression for F then gives simply 

F — Ape + \Bp%, (26) 

where A = A c — \aA C o and \B = B c + a 2 B C o + a 4 B - The transition temperature is 
defined by A = a(T - T c ), so that 

aT c = A + + {A 2 _+AG 2 13 }\ (27) 

where A± = A u + A 12 ± A u and A C o = 4Gi3 h &s been used. 

For later purposes, we note here that minimization of F above yields p 2 c = —A/B = 
(a/B)(T c — T) for the primary order parameter, and that the specific-heat anomaly at T c is 
given by 

AC = C(T C + ) - C(T-) = -t£. (28) 



III. ELASTIC COUPLING 

The orthorhombic symmetry of CuGeO^ allows for nine independent elastic constants: 
Cn, C12, C13, C22, C23, C33, C44, C55, Cm- The elastic energy is conveniently written as 

F* = \Y, G wW ( 29 ) 

where Voigt notation (1 = xx, 2 = yy, 3 = zz, 4 = yz, 5 = zx, 6 = xy) has been used. 

Lowest order coupling between the strain tensor and displacement vector p is of the form 
~ ep 2 and also has nine independent terms 

F Pe = ^ [D^XXPXPX + D 2 e y yPyPy + D 3 C ZZ P Z P Z + D^xxPyPy + tyyPxPx) 

+D 5 {e xx p z p z + e zz p x p x ) + D 6 (e yy p z p z + e zz p y p y ) (30) 
+D 7 e yz p y p z + D 8 e xz p x p z + D 9 e xy p x p y ], 



where the short-hand notation A = A(t) and p a pp = p a (R)p/3(R') has been used. The 
same near-neighbor interactions as considered above for F 2 and F 4fe yield 

&xxPx + D 2 e y yp y + D 3 e zz p z + D±{e xx p y + e yy p x ) 
+D bC e xx p 2 z + D 50 e zz p 2 x + D 6C e yy p 2 z + D eo e zz p 2 y + D 9 e xy p x p y ]. (31) 

Recall that Cu displacements are along the z axis and that O displacements are in the xy 
plane. The interaction constants Di, D 2 , -D4, -D50, Dqo, and D$ thus involve — couplings 
whereas D 3 ,D 5C , and D 6C involve Cu — Cu couplings. Cu — O interactions are zero for 
the same reasons as described in section 2A-5, although those inlvolving density-density 
couplings could be added if desired. With the assumption p x = p y = Pol as well as the 
result from Eq. (24) 

Po = XP 2 c, X ^ « 2 + 2app 2 c , (32) 
the elastic coupling term reduces to 

\K x e xx + K 2 e yy + K 3 e zz + K 6 e xy ]p 2 c 
= H K v e vPci (33) 

where 

K x = (Di + Di)x + 2D 5C 
K 2 = (D 2 + D A ) X + 2D 6C 
K 3 = 2D 3 (D 50 + D 60 ) X 

K 6 = D 9X . (34) 

Note that x is temperature dependent through p 2 c and that it is proportional to the square 
of the density-density coupling constant, x ~ G\ 3 . It is also convenient to write 

K„ = K^ + k^pl, (35) 

where K^ is independent of temperature and depends only on Cu—Cu interactions, whereas 
the coefficient of the temperature-dependent contribution depends only 011O-O inter- 
actions. 
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A. Thermal expansion and elastic constants 

The full free energy now takes the simple form0 

F = Ap 2 c + \Bp% + \ C ^e v + £ K^pl. (36) 

flU fl 

Minimization with respect to the strain tensor yields the following results for the effects of 
Cu and O displacements on thermal expansion, the expansion coefficient, and the elastic 
constants 

2 1 „2 



e M = - s^[K vo + k u p 2 c ]p 



c 



Aa„ = -(a/B) J2 s»v[Kuo + 2Kp 2 c ] 
AC, U = (K vo + K&Wvo + KpI)/ b (37) 

where is the compliance matrix (s^ = C^), Aa^ = a M (T+) — « M (T~), and with a 
similar definition for AC^. Due to the temperature dependence of p 2 Cy thermal expansion is 
not simply proportional to the square of the order parameter, and the expansion coefficient 
does not exhibit a simple step discontinuity, as is usually the case with low order (~ ep 2 ) 
interactions. In addition, the elastic constants exhibit temperature dependence below T c , 
also due to the displacement of O ions (since k^p 2 c ~ p 2 ). Such effects can also be mimicked 
by including higher-order couplings. 



B. Numerical fits 

With available experimental data, it is possible to fit a number of the Landau model 
parameters. From the specific heat data of Ref. [IT), the estimate AC ~ —4.5 x 10 4 J / (Km 3 ) 
can be made. Comparing this result with Eq. (28), and using T c = 1A.2K, gives a 2 / B ~ 
3.2 x 10 3 J/m 3 (where the free energy F is in units of energy /volume). The order parameter 
at low temperature takes the form pc(0) = (aT c /B)^. This expression can then be used 
with the above result, and the estimate from Hirota et al. of pc(0) — 0.0014c, to give 
a ~ 2.6 x 10 29 J /(Km 5 ) and B ~ 2.2 x 10 55 J /(Km 7 ) in (cumbersome) SI units. 
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There are also available limited data on the elastic constants near T c . In particular, from 



the results of Ref. |14] on Cn, C22, and C33, it is clear that only C33 exhibits a detectable 
discontinuity so that K w ~ 0, if 20 — 0. From these data, the estimate AC33/C33 ~ 2.5x 10~ 3 
can be made. Using this result, along with@ C33 ~ 3.7 x \Q ll N/m 2 , as well as Eq. (37) 
and the above value for B, the estimate K 30 ~ 1.5 x 10 32 J/m 5 can be made. These results, 
together with the relations (37), imply that the only other elastic constants which could 
display a discontinuity at T c are C 3 q and Cqq, provided that K 60 has appreciable magnitude. 

The temperature dependence of AC^ below T c from Eq. (37) is of the form ~ ACo + 
c iPc + c "iPci where in general pc ~ (T c — Tf. From the data of Ref. |14], there appears to 
be a region near T c where ACn, AC22, and AC33 increase nearly linearly with decreasing 
temperature, followed by a less rapid increase at lower temperatures. A more detailed 
analysis of these results^! suggests a value j3 = 0.42(2). This estimate may be compared 
with p = 0.25(5) found in an earlier study of elastic properties.0 A discussion of possible 
critical behaviors is deferred to the next section. 

It is also of interest to examine the expression for the effective fourth-order coefficient 
B 1 = B—Ys^v K^qKvq which results from elastic coupling.0 With the reasonable assumptions 
that the dominant contribution is from i^o, and that S33 ~ I/C33, one finds (B' — B)/B ~ 
10~ 3 . The argument given in Ref. [l^ that elastic coupling could be responsible for reducing 



the magnitude of the fourth-order coefficient to near zero appears to be inconsistent with 



the above analysis. Thus, the proposal in Ref. [T(] of tricriticality is not supported by the 
present work. 

Experimental results on the three diagonal components of the strain tensor may be found 
in Ref. [17|. From Eq. (37) and the approximations K w = K 2 q = 0, temperature dependence 
of the form 

e M = -[s^K 3 o + (sniki + Sfakz) + s^ 3 k 3 )p 2 c ]p 2 c , (38) 

is expected. It is clear that such predicted behavior is consistent with the data. Detailed 
fitting to extract estimates of the three parameters fc M must await determination of the 
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absolute values of the other elastic constants (only Cn, C22, and C33 are presently available) 
and hence knowledge of the compliance matrix. With the assumption of quadratic order- 



parameter dependence, the authors of Ref. [H] extracted estimates 2/3 = 0.61(2) — 0.69(8). 

Discontinuities in Aa^ at T c along the three principle axes a, b, and c were found to be 
1.3(1), -2.3(1), and -0.4(1) x 10~ 5 /K, respectively^ From (37), as well as the approx- 
imation S33 ~ I/C33, the c-axis expansivity is given by Aa c ~ —{a/B)K^/C^- Using 
the parameter estimates given above, this yields the model prediction Aa c ~ —5 x 10~ 6 , 
in very good agreement with the data. Comparison between the model predictions 
Aa a = — (a/ B)KsoS\3 and = —(a/B) -K30S23 at T c with corresponding experimental 
results requires knowledge of the compliance matrix. It is also of interest to note that the 
temperature dependences observed below T c in the expansivity are accounted for by Eq. 
(37). 



IV. CONCLUSIONS 

In spite of the great interest in the novel properties of the new inorganic spin-Peierls sys- 
tem CuGeO% as reflected by the large number of recently published experimental studies, 
there has been relatively few theoretical investigations. This is partly due to the complex 
nature of the relevant interactions caused by the non-negligible interchain interactions in 
this not-so quasi-one-dimensional system. Phenomenological models are thus particularly 
useful due to their simplicity and serve to guide further study of a more microscopic nature. 
The Landau-type free energy studied in the present work was constructed from symmetry 
arguments with dominant close-neighbor interactions included. The analysis demonstrates 
that the simplest of interactions can account for the model of rather complicated ionic dis- 
placements proposed in Ref. |8]. O-ion distortions are coupled linearly to the principal Cw-ion 
displacements by charge-density interactions. Elastic distortions as a function of tempera- 
ture below T c , not present in simple Landau models, are seen here to be a consequence of 
this incipient O-ion movement. This contrasts with the elastic anomalies at T c , which are 
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due to the Cu displacements. 

A simple argument based on the results of this work suggest that the transition should 
belong to the standard XY universality class (where, e.g., a = —0.012(10), (3 = 0.349(4), v = 
0.671(5)). One component is simply the Cu displacement along the c axis, p z = pc, and 
the other is the concomitant O distortion in the ab plane p± = po- Note that there is not 
rotational degeneracy in the ab plane due to the second-order anisotropy terms: A specific 
direction is favored. 

In addition to the wide ranging estimates for the critical exponent (3 = 0.25 — 0.42 
determined from the elastic properties discussed above, there have been disparate results 
from other techniques. Analysis of neutron diffraction data for (3 and v indicated a simple 
mean-field transition at T C J1! in agreement with one recent study of the specific heat.il In 
contrast, an earlier specific-heat measurement yielded an estimate a ~ 0.4,0 whereas the 
most recently published result claims —0.15 < a < 0. 15.il The above results for j3 may 
also be compared with those from X-ray scattering data:0@ (3 = 0.26(3) and (3 = 0.37(3). 
One may conclude from all of these results that the most recent of studies of the critical 
behavior appear to favor a scenario where the spin-Peierls transition in CuGeO^ belongs to 
one of the standard universality classes (Ising, XY, or Heisenberg). It is of interest to note 
that some of the above results suggest tricitical behavior (a — = |). Simulations of 
simple microscopic models of structural phase transitions indeed indicate that the fourth- 
order coefficient of a corresponding Landau-type expansion can be rather small.il Analysis of 
recent X-ray data on CuGeO% demonstrates that the critical region is rather small and that 
tricritical behavior is evident if exponents are estimated using larger values of the reduced 
temperature.il Similar conclusions were recently made based on Monte Carlo simulations of 
a magnetic system that was previously thought to be tricritical.@ 
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FIGURES 

FIG. 1. Schematic showing displacements of Cu (open circles) and O (filled circles) ions in the 
ab (top planel) and cb (bottom planel) planes below T c , after Ref. 8. Not shown are the two O 
and two Ge ions per high-T unit cell (rectangles) which are undisplaced. 
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